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We discuss the leading order correction to the equation of motion of a particle with spin on an 
arbitrary spacetime. A particle traveling in a curved spacetime is known to trace a geodesic of the 
background spacetime if the mass of the particle is negligibly small. For a spinning particle, it is 
known that there appears a term due to the coupling of the spin and the Riemann tensor of the 
background spacetime. Recently we have found the equation of motion of a non-spinning particle 
which includes the effect of gravitational radiation reaction. This paper is devoted to discussion 
of a consistent derivation of the equation of motion which is corrected both by the spin-Riemann 
\ coupling and the gravitational radiation reaction. 

I/-) . PACS number(s): 04.30.Db, 04.25.Nx 
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G\ ■ I. INTRODUCTION 

The gravitational waves from an inspiralling compact binary is one of the most promising sources expected to be 
detected by the near- future interferometric gravitational wave detectors such as LIGO, VIRGO, GEO, TAMA and 
£jfV LISA In order to extract the information of a binary from the last inspiralling stage, we need accurate theoretical 
templates of the gravitational wave forms [gj. As an approach, perturbations of a black hole by an orbiting point 
particle have been studied 

Most of the previous papers considered the case that the orbiting point particle is structureless and is characterized 
only by a mass parameter. In those papers, the point particle is denned by using the Dirac delta function, therefore 
the induced perturbation diverges at the particle. This leads to some conceptual problems in the perturbation study; 
1) how far we can incorporate the strong non-linearity around the particle, and 2) how definitely we can define the 
motion of the particle though the perturbed metric diverges at the particle. 

In a previous paper [5| we have considered these problems in the case of a non-spinning particle, extensively 
using the technique of the matched asymptotic expansion that had been developed by many authors (e.g., D'Eath 
Q and Thorne and Hartle [Q). The matched asymptotic expansion is a technique to match physical quantities 
in different perturbation schemes. In the case of a non-spinning particle, the procedure goes as follows. We first 
construct the metric in two different regions. We assume the metric around a particle (body's neighborhood 0) to 
be approximated by a Schwarzshild geometry plus its perturbation describing tidal distortions due to the background 
curvature (internal scheme |5j or body expansion §). In this region the metric is expanded in powers of r/L, where r 
is the circumference radius measured in the local inertial frame of the background and L is the characteristic curvature 
scale of the background spacetime. The metric far from the point particle (external universe (?]]) is approximated by 
a given background spacetime plus the perturbation generated by a point source (external scheme |3| or external- 
universe expansion |7]]). The metric in this region is expanded in powers of Gm/r, where m is the mass of the 
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particle. The particle is presumed to be compact and isolated in the sense that Gm is far smaller than L. Then the 
metrics constructed in different expansion schemes are matched in the overlapping region (buffer region JtJ) where 
both expansion schemes are valid. In practice, the metrics are expanded in powers of both r/L and Gm/r and they 
are matched by imposing appropriate coordinate conditions. As we take a Schwarzshild geometry as the lowest order 
approximation of the solution around the particle, non-linearity can be fully taken into account there. 

In ||, the equation of motion has been derived from the condition that the above asymptotic matching of the 
metric can be consistently done. In the internal scheme, the background spacetime is a Schwarzschild geometry which 
has spherical symmetry. This enabled us to decompose the metric perturbation in the internal scheme in terms of 
the tensor spherical harmonics. The important step was to require that the translational gauge modes should vanish 
in the internal scheme, which corresponds to fixing the center of mass of the particle. Equating each component of 
the metrics derived in the internal and external schemes in the overlapping region, we have obtained the equation 
of motion of the center of mass in the backgound spacetime in the external scheme which includes the effect of 
gravitational radiation reaction. It should be noted that Quinn and Wald have recently obtained the same equation of 
motion for a non-spinning particle by an axiomatic approach without resorting to the matched asymptotic expansion 

P. However, it is beyond the scope of this paper to discuss the physical meaning of their axioms in the framework 
the matched asymptotic expansion. 

Having shown that the asymptotic matching can be consistently done with the external metric generated by a non- 
spinning point particle, we have clarified that the expression of the energy momentum tensor for a point particle used 
in computing the gravitational radiation reaction to the orbit in previous works Q can be regarded as a Schwarzschild 
black hole or any compact object provided that it is kept sufficiently spherically symmetric. Since it is unnecessary 
to evaluate the external metric at the location of the particle, we do not encounter the divergence even if we use 
the delta function source. Applying the asymptotic matching, we obtain the equation of motion of the origin of the 
internal frame with respect to the background spacetime. And hence the equation of motion is understood in a well- 
defined manner. As a result, we found that the strong equivalence principle holds to O ((Gm/I) 2 ) for a non-spinning 
compact body in the limit Gm/L <C 1. Here compact body means that it has no typical length scale other than its 
Schwarzshild radius. 

On the other hand, the post-Newtonian calculation shows that the spin of the particle modifies the gravitational 
wave forms of a realistic binary system M, which motivates us to study the problem of incorporating a spinning 
particle into the perturbation theory. For this purpose, we need to know the correction to the equation of motion due 
to spin and the energy momentum tensor of a spinning particle. 

The effect of the spin to the equation of motion was discussed by many authors in various ways j[(J[ll) . Dixon |Tl| ] 
found the covariant description of the energy momentum tensor of an extended body and the equation of motion of 
it. The resultant equation of motion shows that the coupling of the spin and the Riemann curvature affects the orbit 
of the body. However, the Dixon's equation of motion is written in terms of the full metric which can be obtained 
only after the Einstein equations are completely solved. On the other hand, what we are interested in is to describe 
the motion of a compact body and the effect of radiation reaction in a given external background. For this purpose, 
further reduction is necessary. Furthermore, Dixson's equation of motion does not seem to be applicable to a strongly 
self-gravitating body like a black hole in a strict sense. 

The effect of the spin was also discussed, for example, in D'Eath |(| and Thorne and Hartle jjj, in which the strong 
self-gravity of the body is taken into account by using the technique of matched asymptotic expansion. Especially 
Thorne and Hartle |?J established the framework of matched asymptotic expansion of the metric around the particle. 
Using the metric thus derived, they integrated non-covariant conservation laws written in the inertial frame of the 
background spacetime (body's local asymptotic rest frame in on the world tube enclosing the particle. They 
found that these integrations contain the 'laws of motion and precession' which can be converted into the equations 
of motion and spin when an approximate solution of the Einstein equations of the surrounding spacetime is given. 
The laws of motion indicate that the geodesic motion is corrected by the coupling of the spin of the body and the 
Riemann curvature of the surrounding spacetime. As a specific example of a system, they considered the case of 
two Kerr black holes orbiting each other. They assumed that the separation of these two black holes is large enough 
and derived the equations of motion and spin which are accurate up through the 1.5th post-Newtonian order. Since 
the radiation reaction force is known to give rise from 2.5 post-Newtonian order, the correction due to the radiation 
reaction derived in |5| did not appear in [ ?J . 

The purpose of this paper is to derive the correction to the equation of motion due to the radiation reaction as well 
as to the spin-Riemann coupling in a unified manner. Again we use the technique of matched asymptotic expansion 
in order to construct the metric. Recently we have discussed the perturbation due to a spinning point particle |l2| | 
and calculated the gravitational waves from a spinning particle orbiting a Kerr black hole up through the 2.5th post- 
Newtonian order. In this calculation, we have assumed that the spin contribution to the energy momentum tensor is 
given by a derivative of the delta function. The result agrees with the one obtained by the standard post-Newtonian 
calculation J9| in a suitable limit. Hence, we assume that metric far from a spinning particle is also approximated by 
a linear perturbation of a general background spacetime generated by the delta function source, and we assume that 
the metric around the particle is approximated by that of the Kerr geometry with perturbations, which should/ will 
be justified by the consistency condition of the matching. 

The next step is to extract the information of the motion of the particle from the metric thus constructed. At 
this point we have two methods to obtain the equation of motion; the consistency condition of the matching with an 
appropriate choice of gauge in the internal scheme [^J, or the use of non-covariant conservation laws in the inertial 
frame of the background spacetime H. The former method developed in p] seems difficult to apply to the present 
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case. This is because we do not know a method to distinguish the translational gauge modes of the metric perturbation 
in a Kerr background. Hence we cannot fix the center of mass condition in the internal scheme when matching the 
metrics in the internal and external schemes in the overlapping region. Thus we adopt the latter method and integrate 
the conservation laws as in JjJ to derive the equation of motion. Relating to the difficulty of the former method, we 
cannot verify the consistency of matching to the same extent as we could in the case of non-spinning particle. Thus 
we simply assume it. 

This paper is organized as follows. In section 2, we construct the metric in the overlapping region (buffer region), 
using the covariant expansion method of the tensor Green function H. We suppose the reader is familiar with the 
concept of 'bi-tensors' |l3| and skip the computation of the tensor Green function 0. In section 3, we derive the 
equation of motion of the spinning point particle. Section 4 summarizes the result. 



II. MATCHED ASYMPTOTIC EXPANSION 

Following the procedure given in || , we first discuss the external scheme of the metric. We have a general background 
metric g^ v {x) which satisfies the vacuum Einstein equations around the particle in the external scheme. We compute 
its linear perturbation with the Green function metho d. T he perturbation in the external scheme is constructed with 
an assumption on the perturbation source. (See Eq. (2.3) below.) We will discuss the consistency of this assumption 
later when we compute the metric in the overlapping region. 

We put the metric gy. v {x) + h^ v {x) into the Einstein equations and solve the linearized equations for h^ v assuming 
that the background metric <7 M „ satisfies the vacuum Einstein equations around the particle of interest^. We introduce 
the trace- reversed metric perturbation ^> M „(x), and put the covariant harmonic gauge condition on it. 



Vv(x) = h^ v {x) - -g^(x)g^ (x)h^(x) 







(2.1) 
(2.2) 



We can make use of the tensor Green function of the linearized Einstein equations with the retarded boundary 
condition derived in Sec. 2 of 0]. Here we assume the metric perturbation is generated by the following source term 
0, 



T^(x) =m dTl v»(x, T)v v {x, T) 



S^(x-z(T)) 



S^(x,T)v v \x,T) 



5W(x-z(T)) 



(2.3) 



where v fl (x,T) = g^ a {x,z{T))z a (T) and S^{x,T) = g» a (x, z(T))g u (x, z{T))S al3 (T). g» a (x,z(T)) is a bivector of 
parallel displacement defined, for example, in Appendix A of j5j, and S a ^(T) is an anti-symmetric tensor which is 
called the spin tensor of the particle and is assumed to satisfy the center of mass condition, S a0 (T)v^(T) = 0. A 
dot over a function means T-derivative of it, such as z a (T) — dz a / dT. The amplitude of the spin tensor S a p(T) is 
usually assumed to be of order m in the case of a Kerr black hole, which seems to indicate the 2nd order perturbation 
theory is necessary when we want to discuss the spin effect of the particle in a consistent expansion with respect to m. 
However, we will find it not so if we assume that the construction of the metric by the matched asymptotic expansion 
is consistent, terms of order m 2 , the 2nd order perturbation theory is unnecessary. The exact form of the metric 
perturbation becomes 



^ v (x) = 2Gm\ 



1 



^(x,z(T))z a (T)z?(T) 



&(x,z(T)) 



rT{ , : { r) \u^ a0 (x,z(T))a n (x, z(T))S-> a (T) Z P(T) 



a 3 (x,z(T)) 



1 



&(x,z{T)) 

1 d 



a 2 (x,z(T))dT 



u^ a ^(x,z(T))S^ a (T)z (T) 

{u^ aP {x, z(T))a. n (x, z(T))S~< a (T)zJ(T)) 



§In this derivation, we follow the notation used in jl2 
z(T), and the indices fj,,u,£,p for the field point x 



We assign the indices a, (3, 7, 8 for the point on the particle trajectory 
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&(x,z(T)) 

T Ret (x) 



v» v a p(x, z(T))a ;i (x, z{T))S^{T)z {T) 



T=T Rct (x) 

dT (v^ af3 (x, z(T))z a (T) Z P(T) + v^ af3;l (x, z(T))S^ a {T)z' 3 {T))^ , (2.4) 

where cr(x, z) is a half the squared geodetic interval between x and z and is defined, for example, in Appendix A of 
u MI/a/3 (x, z) and v tJ "" a P(x, z) appear in the expression for the tensor Green function and are defined in Sec. 2 of 
and Tn et (x) is the retarded time of the particle which is a scalar function determined by a (x, z(Tji et )) = with 
the condition that z{T^ et ) is on the past light cone of x. 

We next consider the metric in the overlapping region, in which the metric is expanded in powers of Gm/r and 
r/L. We first introduce coordinates {T, X 1 } = {X a }, (i = 1,2, 3; a = 0,1,2,3)**, such that they become those 
introduced in fl3]| when taking the test particle limit m — > : 

lim (a.. a (x,z(T)) + e m (T)X l ) = 0, (2.5) 
limi 2 (T) = -l, (2.6) 

m— i-0 

lim z a {T)e m {T) = 0, (2.7) 

m— *0 

lim e a l {T)e a3 {T) = (2.8) 

m— >0 

We assume that X % X % jl? is small enough in the overlapping region so that the coordinate transformation x a = x a (X a ) 
is single valued Eqs. ( |2.7| ) and ( |2.8| ) show that e a i(T) are the spatial triad at z(T) on the 3-hypersurface normal 
to i(T). The coordinate system for finite m shall be determined later. 

We denote the components of the metric whose order of magnitude are 0({Gm) n / L m ) by ^\ /i a &. 

oc 

JS= E ( 2 - 9 ) 

m,n— 

We call these terms the (™)-components of the metric Simple dimensional analysis shows 

{ ™ } h ab ~\X\( m - n \ (2.10) 



where = yX^X 1 . In the next section, we use the Landau-Lifshitz pseudotensor in order to compute the equation 
of motion, and the metric is expanded by taking the flat Minkowski metric r\ a b as a background. For its purpose, we 
raise the indices of the (™)-components of the metric by rj a i,, and we define the trace-reversed (™)-components of the 
metric with respect to the flat Minkowski. 

IJ^g^^wl?^ (2 ' n) 

The (™)-components of the metric are obtained from the background metric g^^x). The background metric in the 
coordinates {X a } is partly given in H. A further calculation results in 

g„ v (x)dx»dx» = (V(T) + 2z"(T)^i(T)X i + ^*(T)^(T)X*^ 

+R a01 s(z(T))z a (T)XP(T)z-<(T)X s (T) + 0(|X| 3 ))dT 2 

+2(V(T)e m (T) + e a i{T)^-{T)Xi 

+ ^R af 3 1 s(^T))e a l (T)Xf } (T)P(T)X s (T) + 0(\X\ 3 )) dTdX* 



** We assign the indices, a,b,c,d, for the spacetime coordinates i.e. 0,1,2,3, while the indices, i,j,k,l, for the spatial 
coordinates i.e. 1,2,3. Minkowski and Kronecker summention conventions are taken. 
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+ (e Q l (T)e cy (T) 

^i?^ 75 (z(T))e Q l (T)X' 3 (T) e ^(r)X 5 (r) + 0{\X\ 3 )^jdX l dX^ 



where X a {T) = e a t (T)X l . Because of ( |2.10| ), we find 

° ]h 00 = z 2 (T) + 0(Gm/L), 
§hoi = i a (T)e m (T) + 0(Gm/L), 
l]h. tJ = e a l {T)e a] {T) + 0(Gm/L), 



Dp 

2z a (T)-^(T)X* + 0(Gm\X\/L 2 ), 
l]h Ql = e a t (T)^(T)XJ + 0{Gm\X\/L% 



dT 



^hij = 0(Gm\X\/L 2 ) } 



oj^oo = R a pjs(z(T))z a (T)X 13 (T)z' 7 (T)X 5 (T) + 0(Gm\X\ 2 / L 3 ), 
oj/KH = ^ Q /3^(2(T))e Q l (r)X' 3 (T)iT(T)X 5 (T) + 0(Gm\X\ 2 / L 3 ), 



fli'i — 



i J R Q/37 5(z(T))e Q l (T)X' 3 (T)e^(r)X 5 (r) + 0(Gm\X\ 2 / L 3 ). 



(2.12) 

(2.13) 
(2.14) 
(2.15) 

(2.16) 

(2.17) 
(2.18) 
(2.19) 

(2.20) 
(2.21) 



As was argued in 0, the matching condition requires that the (g)-components of the metric ( |2.13 )-( 2.15) should 
be equal to those of the flat metric 



(0) 

(0) 



hab = Vab 



(2.22) 



which is exactly realized by the coordinate conditions ( 2.6 )-( 2.8 ) . 

The ((^-components of the metric should vanish in order to put the particle in the locally rest frame of the 
background geometry. 



From (2.16)-(2.18), we thus obtain 



(0)" ab 



z a (T)^p(T) = 0(Gm/L 2 ), 

Dp 

e a t (T)-^(T) = 0(Gm/L 2 ). 



From (2.6), (2.7) and ( |2.24 ), we get the geodesic motion of the particle in the background geometry. 

D 



dT 



z a (T) = 0(Gm/L 2 ) 



(2.24) and (2.25) show the geodesic transform of the spatial triad e a i(T) in the background geometry. 

De a t 



dT 



■(T) = 0(Gm/L 2 ) 



(2.23) 

(2.24) 
(2.25) 

(2.26) 

(2.27) 



The (™) -components of the metric and the (™)-component s of the metric due to the spin of the particle come from 
both the background metric g M „ and the metric perturbation (12.41 ). The covariant expansion of the metric perturbation 
due to the monopole particle is already given in [0. Using fl2.5| ), we can soon get the expression of it up to 0(mX). 
We further compute the metric perturbation due to the spin up to 0(mS\X\°). 
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''i- { - r) ' Lr '' ' Lr " = i 2(; '"'( m + ^\^'^ T )) ia ( T ) x0 ^ i ^ T ) xS i T ) 



2 {V a p(T) + V aPl (T)X\T)) z«(T)zP(T) + 0(m\X\ 2 , m\S\\X\") SdT 2 



-2l -2Gmi? a/ 3 75 (z(r))z Q (r)X' 3 (T)i''(r)e 5 l (T) 



-2 (V afj (T) + V a01 {T)X~<(T)) z a (T)e^(T) 



-2Gj^S aP (T)X a (T)e^(T) + O( m \X\ 2 ,m\S\\X\ ) { jdTdX t 
+ \^-AGm\X\R aPlS {z(T))z a {T)e^{T)P{T)e s J {T) 

-2 (V aP (T) + V aM {T)X^{T)) e a t (T)e^(T) + 0(m\X\ 2 , m\S\\X\°) \dX l dX^ 



V af3 (T) = Gm / dT'(v a/3 a , ,(z(T),z(T'))z a '(T')zP'(T') 

J —OG 

+« a/ W ; y (z(T), z(T'))S-'' a ' (T')zP' (TO) , 
V a ^{T) = Gm f T dT'fv^a'pt.^ziT)^^'))^' (T')z 0> (T') 

J —OO 

+^ /3 a' / 3'; 77 '(^(T),z(T'))^' Q '(T')i /3 '(T')). 



We first discuss the (^-components and the (^-components of the metric. From ( 2.28 ), we obtain 



(0)7 
(!)"00 


AGm 

\x\ ' 


(0)7 






= 0, 


(0)7 
(2) "00 


spin 


-o, 


"(0)7 
(2) /l 0i 


spin 


2Gm 

\X\3 




spin 


= 0. 



S af3 (T)X a (T)eMT) + 0((Gm) 3 /L\X\ 2 ), 



(2.28) 

(2.29) 

(2.30) 

(2.31) 

(2.32) 
(2.33) 
(2.34) 

(2.35) 
(2.36) 



where [ ] sp i n means that only the spin contribution is taken into account. In view of the internal scheme, the (°J- 
components of the metric describe the background metric, thus the consistency condition of the matched asymptotic 
expansion requires that these components realize the asymptotic gravitational field of a Kerr black hole. |g| argues 
the (5)-components of the metric fit the asymptotic metric of a Schwarzschild black hole with mass parameter m and 
we have shown in |L2f that the spin contribution to the (^-components of the metric correctly realizes the asymptotic 
field of a Kerr blacK hole. Thexefpre the metric thus constructed in the external scheme matches the metric in the 
internal scheme and the use of ( |2.3[ ) is justified in computing the leading order contribution of the mass m and the spin 

S. The monopole contribution to the (^-components of the metric cannot be computed in the external scheme unless 
we develop the 2nd order perturbation formalism. However, the consistency condition of the matched asymptotic 
expansion requires that these components are those of the asymptotic expansion of a Schwarzschild black hole. Thus 
we finally obtain 



(0) 



hnn — 



(2) "00 



(0) 



hni — 



(2)"0i 



(Gm) 2 
\X\* 
2Gm 
W 



S a0 (T)X a (T)e^(T) + 0((Gm) 3 /L\X\ 2 ), 



(2.37) 
(2.38) 



G 



(0)-, _ ( G ™) 2 ( O* , X ' X] 



(^-^(-26 i] + lw j. (2.39) 

We next discuss the (})-components and the (^-components of the metric. As is argued in 0, we can take the 
coordinate system such that all the (^)-components vanish as long as the (^-components of the metric vanish. In view 
of the internal scheme, the (^)-components of the metric describe a linear perturbation of the background metric. 
This perturbation is induced by the curvature of the external universe, i.e. by the (J)-components. However, we 
have seen that the (J)-components should vanish because of the coordinate condition. Thus the (^)-components of 
the metric have no physical mode, and we can set them zero with a suitable choice of coordinates. As discussed by 
Thorne and Hartle 0, it is essential to put them zero for reducing the uncertainties in the definitions of the mass, 

momentum and spin of the particle when deriving the equation of motion. Hence we require /„\/i a 6 — 0. From (2.12) 

and (2.28), we obtain (^-components of the metric. 

S^oo = \ (e a i(T)e ai (T) + z 2 (T) - 2) 

-2V a/} (T)z a (T)z (T) + 0{{Gmf/L 2 ) (2.40) 
g/iw = z a (T)e ai (T) 

-2V aP (T)z a (T)e^(T) + 0((Gm) 2 / L 2 ) (2.41) 
gjSy = {e a l {T)e aj {T) - <%) - i<% (e a k (T)e ak (T) - z 2 (T) - 4) 

-2^(T) e Q 4 (T)e^(T) + 0((Gm) 2 /L 2 ) (2.42) 



We then obtain the coordinate conditions which satisfy |Jjft. a fc = 0. 

( 7 ;a (x,z(T)) + e ai (T)X i = (2.43) 

i 2 (T) = -1 + V aP {T) (g al3 (z(T)) + 2z a (T)zP(T)) + 0((Gm) 2 /L 2 ) (2.44) 

z a (T)e m (T) = -2V al3 (T)z a {T)e f \(T) + 0((Gm) 2 /L 2 ) (2.45) 

e a i(T)e aj (T) = % + V a p{T) (S l3 g al3 (z(T)) + 2e tt l (T)e^(T)) + 0((Gm) 2 /L 2 ) (2.46) 

Setting the (^-components of the metric zero will give us the coordinate conditions to the order 0((Gm) 2 / L 2 ), but 
we do not need them for the present purpose. 

The metric derived by the matched asymptotic expansion is summarized as follows. 

[oj/laft = Vab , (2.47) 

gji«b = 0, (2.48) 

g/H.0 = 2 -R a p, S (z{T))z%T)xP(T)P{T)X\T) 

+0(Gm\X\ 2 /L 3 ), (2.49) 

$h oi = ^ Q/375 (z(T))e Q l (T)A' 3 (T)^(T)A 5 (T) 

+0(Gm\X\ 2 /L 3 ), (2.50) 
1 

+0(Gm\X\ 2 /L 3 ) , (2.51) 

(0) 7 _ 4GTO (9 

(1) 00 ~~ ~~\X\~ ' (<Z-OZ) 
P ) h oi = 0, (2.53) 
[°]^ = 0, (2.54) 
( $h ab = 0, (2.55) 
[l]^oo — z a (T)-^e a i(T)X i (T) 



S% = -R a ^s(z(T))Xf'(T)X s (T) (e a t (T)e^(T) + S l3 z a (T)r(T)) 
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+ ^^R a ^s(z(T))z«(T)XP(T)z~<(T)X s (T) 

-2V a[il {T)i a {T)iP(T)X~t{T) + 0((Gm) 2 \X\/L 3 ) , (2.56) 
gftoi = e a i(T)-^e aj (T)Xi 

+2GmR aPlS (z(T))XP(T)P(T) (- z a (T)e 6 \(T) + pe«,(T)I s (T)J 
-2^ 7 (T)i Q (T)e' 3 4 (T)^(r) + 0((Gm) 2 \X\/L 3 ) , (2.57) 
W- hij=StJ z°(T)^e ak (T)X k 

+2GmR afJlS (z(T))[^X^(T)X s (T) (e a i(T)e l j(T) — d i3 z a (T)z 7 (T)) 
-2\X\z a (T)e fj l (T)z'-'(T)e s :j (Tf 



-2V a ^(T)e a i(T)e^(T)X^(T) + 0((Gmf\X\/ L 3 ) , (2.58) 
(0)7 _ (Gm) 2 

(2) /l oo - | X | 2 , (^-oyj 

$hi = -^S a0 (T)X a (T)e^(T) + 0((Gm) 3 /L\X\ 2 ) , (2.60) 



(^^wr^wJ ' (2 - 61) 

[2)^6 = 0. (2.62) 



III. EQUATION OF MOTION 

We apply Eqs. (2.2a), (2.2b), (2.3a) and (2.3b) of jjj to derive the equation of motion, which are 

1 

16ttG 



po ( T > r ) = 77T75 / d 2 ^tf ab0j , & (3.1) 

\X\=r 



J« (T, r) = — !— / rf 2 ^ (X l W Mk a - x :j H ia0k a + H lk0j - H jm ) (3.2) 
16ttG J\x\=t 

±P*(T,r) = -± J s d 2 S^{X) (3.3) 

J" (T, r) = ~f_ d 2 S k (XH^(X) - XH ik (X)) (3.4) 



s 

As in Q, we adopt the choice of Landau-Lifshitz for H abcd and t ab . 

jjabcd = jjabcd ^ (35) 

^ = (-g)tf_ L . (3.6) 

The explicit expressions are given in standard texts, such as, Sec. 20 of Misner-Thorne- Wheeler Q] or Sec. 96 of Landau- 
Lifshitz p5| . The integrals are taken over a 2-surface of constant \X\ in the overlapping region for definiteness, but 
we focus on the r-independent terms in the integrals since it is enough only to consider them in order to derive the 
leading order correcti on to the orbit z(T). 

Before considering (EOT) and (|3.3|), we first show that the spin ten sor S a 3(T) is geodetic transported in the back- 



ground geometry g a g(z(T)) in the test particle limit m — > 0. Eq. (3.2) has a dimension of (mass) x (length) and 
we consider the terms of order (Gm) 2 . Power counting of X shows that there will be contributions linear in the 
(^-components of the metric and those from bilinear combinations of the (°)- and (^-components of the metric. We 
obtain 



J ij (T,r) = mS a0 (T)e ai e m + 0(G 2 m 3 /L) + (r-dependent terms). (3.7) 
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Eq. (3.4) has a dimension of (mass) 1 and we consider the terms of order Gm 2 /L in the same way. Power counting 
of X shows that there will be contributions from bilinear combinations of the (J)- and (°)-components of the metric, 
and we soon find that (3.4) vanishes. 



— J tj (T,r) = 0(G 2 m 3 /L 2 ) + (r-dependent terms) 
dT 



(3.8) 



Since the spatial triad are geodetic transported in the background geometry to the leading order, fl3.7| ) and ( |3.8| ) result 
in 



DS 



a/3 



dT 



-(T) = 0. 



(3.9) 



We next consider (3.1) and (3.3). Eq. (3.1) has a dimension of (mass) 1 and we consider the terms of order m 
and Gm 2 / L. We find that there will be linear contributions from (5)-, (^-components of the metric, and bilinear 
contributions from pairs of (?)— and (^-components of the metric. We obtain 

P°(T, r) = m + 0(G 2 m 3 /L 2 ) + (r-dependent terms), 



P l (T,r) = 0(G z m 6 IL*) + (r-dependent terms). 



(3.10) 
(3.11) 



Eq. (3.3) has a dimension of (mass) / (length) and we consider the terms of order m/L and Gm 2 /L 2 . There will be 
bilinear contributions from pairs of the (§)— and (^-components and pairs of the (5)— and (^-components of the 
metric. The former pairs give the leading correction to the motion due to the spin of the particle, and the latter pairs 
give the radiation reaction of motion. A straightforward computation results in 



— P u (T,r) = 0(G 2 m 3 /L 3 ) + (r-dependent terms). 
dT 

-^P l (T,r) = |i? Q/375 (z(T)) e " l (T)i' 3 (T)5^(r) 

-^V af3l (T)e\(T) (2z a (T)z^T)+g^(z(T))) 

+mz a (T)—e m (T) + 0(G 2 m 3 /L 3 ) + (r-dependent terms). 
dT 



Taking the T-derivative of ( 2.44 ) and ( 2.45 ), we finally get the equation of motion, 

-^i"(T) = Vf3 7 s(T) (2^(T) g ^(z(T))z s (T) - & \T)P '(T)g ai >(z(T))) 
+ l -R a fi ~ lS (z(T))zP(T)S^(T) + 0(G 2 m 2 /L 3 ), 
where we have defined V a p 7 = V r Q/ 3 7 — l/2g a pV s s y . Introducing the proper time r of the orbit, 



dr 

— = l-V af3 (T)z a (T)zP(T), 



we finally get 



D dz a 
dr dr 



(r) = % lS (r) 



dz a dz 13 dz~< dz 5 
dr dr dr dr 
dz 



dz p 

2^-g^(z(r)) 
dr dr 



dz 5 dz' 3 dz^ aS 



dr dr 



9 ad (z(r)) 



1 rj 7 P 

+ -R« 0jS (z(t))—S^(t) + 0(G 2 m 2 /L 3 ). 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



IV. CONCLUSION 



In this paper, we have derived the equation of motion of a spinning particle on a given background spacetime. We 
have constructed the metric using the technique of matched asymptotic expansion and integrated the conservation 
laws introduced in to extract the information of the equations of motion. One possible drawback of the present 
derivation of the equation of motion could be that the consistency of the matched asymptotic expansion has not been 
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explicitly demonstrated. To do so, we need to extract the translational gauge modes from the metric perturbation in 
the internal scheme and examine if the condition of their disappearance leads to the equation of motion, as done in 
the previous derivation of the equation of motion for a non-spinning particle in |^| . Unfortunately, as mentioned in 
Introduction, since we have no theory of the metric perturbation in the Kerr background, it is at present impossible 
for us to perform this procedure. However, since we find no logical flaw in the method of using the conservation laws, 
we believe the present derivation of the equation of motion is perfectly legitimate. 

The conclusion of our previous paper [p] is that the non-spinning particle moves along a geodesic of the regularized 
perturbed space, 

9ixv(x) = g^{x) + h {v)liu (x), (4.1) 



where h( v \ /xl ,(x) is the tail part of the metric perturbation (2.4) 



h( v )nv(%) = ip( v )nA x ) ~ \g^{ x )9 ip {x)^(v)i P {x), (4.2) 

= -2Gm [ TRCt(X) dT (v^ aP (x, z(T))i a (T)z^(T) + v^ a0n (x, z(T))S^(T)z^ >(T)) . (4.3) 

We have found that the spin of the particle gives rise to an effective force in addition to this, which is the same as 
derived in |10|Jll] fl. Hence in terms of the regularized metric g^ v and renormalizcd proper time f, the equation of 
motion coincides with the one derived in |h],[]T].[7| : 

-|i(f ) = \R a Pl s{z{T))if\r)S^{T) . (4.4) 

As was commented in it is still under investigation to find a method to explicitly compute h^ v ^ u (x) even for 
well-known background geometries, such as a Kerr geometry. 
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